Об определенимйм плотностного контакта в сложнопостроенной среде 

Предложенью аналитическиє модели гравитационного поля и геологической средь для горизонтально-слойстой средьт с 
несколькими плотностньми границами раздела слоев в нижнем полупространстве. Первая модель представляєт собой уравнениє 
силь тяжести, полученное из рассмотрения нормального потенциала силь тяжести в локальной окрестности текущей точки из- 
мерений в двумерном представлений. Для модели средь, в дополнение к известньм в теории интерпретации потенциальньтх полей 
основньм классам контактньх поверхностей (Нумерова, Остромогильского и Страхова), введен новьгй класс Черного контактньтх. 
Приведень образцьт теорем разделения полей для случая нескольких односвязньїх обьемов и нескольких непересекающихся слоев (2 
самьїх популярньх начальньх приближения средьт) в изучаємой модели. 
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Введениєе. Сейчас много внимания уделяется вопросам переинтерпретации материала сьбмок ХХ ст., ввиду 
резкого ограничения обьдмов геологоразведки. Но, ВСЛедствие пОВЬІШ'енНЬХ требований к точности и зффекти- 
вности интерпретации, переработка больших обьдмов геофизических данньчаих требует новьіх подходов к числен- 
ному обоснованию методов, используємьх при обработке соответствующего картографического материала. 

Современнье МеТОДЬ обработки геофизических данньхх развиваются с прицелом на созданиє цифровьгх 
аналитических моделей поля и геологической средьт и вьшсокоточньзх алгоритмов, ориентированньх на приме- 
нение в пакетах ГИС. Особенно зто касаєтся отображения входной информации и результатов ее трансформа- 
ции в сопоставлениий с априорной информацией, представленной графически (визуальная интерпретация). 

Такая , послойная" обработка -- новьій зтап качественной интерпретации, в котором возможности представ- 
ления информации в виде злектронного журнала измерений (Якимчик, 2010) можно использовать для предобра- 
ботки информации в среде специализированньзх пакетов программ. Примерь такой предобработки: вьіщделениє 
регионального фона (минимизация невязки между набором данньх и некоторьшм полиномом) в Мапа; предва- 
рительная фильтрация измерений в Ехсеї; нахождение особьх точек поля (метод обратной интерполяции поли- 
нома Лагранжа) в Маріе; редукции поля с помощью алгоритмов В.Н. Страхова в Руфоп, и т.п. 

Кк интерпретации данньтх потенциальньах полей, адекватной современной геофизической практике, предья- 
вляют следующие требования: 

е  вьшсокая точность и технологичность алгоритмов и программ, их переносимость; 

е адаптация моделей поля к заданию исходньах измерений на коротких профилях; 

е адаптация моделей геолог. средь к ее сложному (нелинейному неоднородному) СоСТтоОяЯниЮ; 

Одно из направлений детерминированного подхода к интерпретации - разработка новьтх аналитических конст- 
рукций геофизического поля (Якимчик, 2009) и геологической средь. Среди последних авторам близки рас- 
смотрения плотностньах моделей типа контакта, одну из которьжх предлагаєм вашему вниманию. 

Определение контакта в известньх моделях средьт (Страхова, Нумерова, Маловичко, см. (Чорний, Дубовен- 
ко, 2002)) традиционно замькаєтся на вьтчисление разности ундуляций теоретической модели средьт от некого 
опорного ее злемента, от которого зти ундуляции не сильно уклоняются. Но при зтом необходимо заранее знать 
геометрию (глубину, форму, ориентацию) опорного злемента (модельного тела либо контакта), как начального 
приближения итераций. От вьібора удачного начального приближения зависит в целом зффективность любого 
из существующих итерационньх алгоритмов (Старостенко, 1978). 

Постановка задачи. Наш способ определения контакта можно обобщить таким образом. В горизонтально- 


слоистой среде с п плотностньми границами раздела слобв ОС; и СС, в нижнем полупространстве С СС, «б», 
баси, С, «С, 0дб,, С з ЕЗЛС, дб'з ОС, представленьт однородньєе тяготеющие массьт извест- 
ной плотности 52) зд, (є) - б, (2) - разность плотностей смежнькх слобв, где 2 є С,, 52) -0 при Є є С. 


Для такого описания средьг полученьг аналитические модели поля и геологической средь. Модель поля -- 
уравнение силь тяжести, полученное из 2) рассмотрения нормального потенциала силь тяжести в окрестности 
текущей точки измерений. Размер окрестности зависит от размеров тяготеющей области и точности аппрокси- 
мации нормальной составляющей потенциала (Дубовенко, 2011). 

Для вьтянутьх вдоль простирания аномалий силь тяжести основное уравнение поля силь тяжести равно 


МАР ПА СОВА НО) 
жін св тт й 


где 2, -2(5) и 2, 20,(2), ас с ср - кривье, ограничивающие снизу и сверху тяготеющую двумерную 


область 5. 

Для модели средьв, в дополнение к известньм в теории интерпретации потенциальньїх полей основньтюм 
классам контактньїх поверхностей (Нумеров, 1930; Остромогильский, 1970; Страхов, 1974), введен новьій класс 
контактньх поверхностей. Зто класс Черного - математическая модель контактной поверхности, описьваємая 


так. В горизонтальной полосе С", ограниченной сверху и снизу прямьіми 2, - й, 2-Й, О« ЛП, «Й, «-юо, 
вьшщделеньгш 2 слоя (плотности б, и б, ). Их разделяет бесконечно гладкая кривая СС , звіздная относительно 


бесконечно удаленной точки, имеет асимптоту 2 - й, п, « П « й,, причем 
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«-нх а», Бі»Ф, «Фуєдіг З 


Если значение П и величина скачка плотности Дд з д, - д, известнь, граница 2 С раздела слобв по вне- 


з І . 
шнему полю восстановима однозначно в классе Черного: СП"? (--со, ср) с 517(-о0, ор). 

Поскольку зта модель отображаєт сложное распределение тяготеющих масе (п попарно непересекающихся 
рудньгх тел или границ раздела слобв, рис. 1), то центральной проблемой является единственность решения об- 
ратной задачи гравиметрии. 





ХУУУКХ У ЗВам ох ох З КькАХ АХ 


Рис. 1. Модель контактной средь в классе Черного. 


Единственность имеет место в достаточно узких подклассах упомянутьхїх классов моделей средьї, т.е., при 
наличий большого обьбма априорной информации об источниках. Если зтого обеспечить нельзя, остабтся упо- 
вать на известньтєе следствия из теорем разделения потенциальньгїх полей. Они позволяют свести решение обрат- 
ной задачи гравиметриий по любому внешнему полю от п масс к решению обратной задачи гравиметрии для от- 
дельньх обеектов, по полям от зтих обьектов. Приведеньт теоремьт разделения полей для нескольких односвя- 
зньїх обьдмов и нескольких непересекающихся слобв (популярньгю8е начальньтюе приближения средьг) в изучаємой 
модели. 

Методика решения. В плоском случає (1) сводится к решению нелинейного интегрального уравнения 


(2- х)? -Н2 
Пий 
з Їто ху ДФ 


по отношению к искомой функции 2 2(2). брннені Р коді рассматривать, как функционал на неком 


0) 





линеале (множество допустимьгтх решений). Оно популярно при восстановлениий потенциала силь тяжести для 
вьтянутьхх по простиранию масс, Задача имеет единственное решение на классе Страхова 51/17 (0) (Старосте- 
нко и др., 1992). 

Разлагая уравнение (3) в ряд по 2 з 2(2), вьтічисляя приращения Де(х, 2) 2 2(х, 29-Л)- 8(х, 2), где (а) - Н(Б)-0, 
можно предложить ряд алгоритмов восстановления контакта по заданному полю. В частности, используем 
модификацию метода Ньютона: 


Д «(8) з са (6) - 2, (5), (4) 
і 22.42) -х'зн? - - 
п - Д 2, (2346, 
Поза паз інібє їн ер вкдз Ї со г РАЄ 


, (2) 
либо модификации метода простьшкх итераций типа Контакт" в двух версиях: 


З (5) 2 2, (2)-- й, (5), п-0,1,2,..., (а) -0, 


























240) ча ху 4 Н? - ме, 
- 5 42 Кі а 5) 
І сон ФО се Ї; -х) од НН | Коган 12) он | 
оре зер зо 
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В алгоритмах бю можно кроме искомого злемента 27 - (2) определить асимптоту Н « соп5і, для чего по- 
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лученьг соответствующиєе аналитическиєе вьражения: 


2 С ПИ С 
сосна село Їє й | о ан (9; 


лк (8) З 2,5), (2) На з Н, ЕТ, п-0,1,2,..., В(2)-0, с -0, ака, 

(8-х) з Най 
во зад) 
Є-х) ЕН 








42 - 209); (7) 


о - тн 425-809) (8) 








дао) 
70 і 1 а 2 слаб з 2 2 а6- - 
асо КОЗУ Ї уні нні ово зарубтвіо 
єоху за 1» 
сь са Р (9) 
«Пе ху нір аа ху нє (5 дає 


Проблемьт вьшчислений. В вьгражениях (4)-(9) существенна проблема вьюбора начального приближения 
контакта. Наши соображения таковь. Ввиду вьтбора модели средьт (рис. 1) имеєм дело с многозкстремальной 
кривой аномалий силь тяжести. Определим на ней найбольший максимум и найменьший минимум аномалий 


з 8(с,) з вир(ежи8(0), 8, з 8(с,) з іНехі 8(х)), и то же для приближения контакта 2, 2: 2(с,) - ар ема с) 
їс, Ф) Іс, Ф) Іс, 
2 «С р зіпБехіг с(2)). Затем задабм нулевоє приближение контакта, каждьюй зкстремум которого соответству- 
; Іс, 4) 


ет однойменному зкстремуму поля силь тяжести в силу того, что б »0 , ОДНИМ ИЗ способов 


о г об) 2-Но дей, Но), ассаь, (10) 
о З 2(5)2-Но ЧКоДя(с), акхссаь, (11) 
где Н,» 0 - задаваємая константой ,асимптота? функции 2, на левом конце промежутка |а, РІ, 


Де(С) « 2(2) - є(а), Де(2, Н,) - аналитически продолженное значение Де (2) на уровне - Ні, Ко - К(Н.) - 
определяемьій весовой козффициент. 
Весовье козффициентьт определяем, но из - оба враження (3) в виде 
зн? 


2 По-зооог ЧИ, «ШО зо ой (2) 
Їт и па б | из 


а 5 а 
гдеизб-с, дизаб, аза-с, Б бос. дЗтот интеграл можно получить в квадратурах - исходньтх для 
определения 2, и 2, при заданнькх Н, 8, и 8,. А приближения 7, значений 2., 27, определяєм по схеме Че- 


бьшева (Березин, Жидков, 1960): 





ат, Де, ПС ке; і (13) 
РА з 


где Ї (2) - линейное преобразование полученньх из (12) квадратур. Процесе (13) при удачном вьшборе с, обла- 


М 0 я 
даєт бьістрой сходимостью (3-4 итерации для получения 2. 7 с точностью 107). 

Вьводь. Теоретически зти конструкции показьвают лучшую сходимость по сравнению с классическими 
примерами (Гравиразведка, 1990), но их практическая зффективность требуєт проверки на подготовленном ма- 


териале полевьжх гравиметрических сьемок. Попьтки обойтись грубьіми приближениями си С; (вьтчисленнь- 


ми не по (13)), вообще говоря, успеха не имели, хотя в ряде случаев приводили к хорошим результатам. 
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